Dark Bell states in tunnel coupled spin qubits 
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We investigate the dynamical purification of maximally entangled electron states by transport 
through coupled quantum dots. Under resonant ac driving and coherent tunneling, even-parity 
Bell states perform Rabi oscillations that decouple from the environment, leading to a dark state. 
The two electrons remain spatially separated, one in each quantum dot. We propose configurations 
where this effect will prove as anti-resonances in transport spectroscopy experiments. 
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Introduction — Quantum coherence allows open multi- 
level systems to form superpositions that uncouple from 
the dissipative dynamics. The system evolves towards 
a stationary pure state. These are named dark states 
from their discovery as non-absorbing resonances in il- 
luminated atom gases Two ground states coupled 
to the same excited state form a superposition that is 
not affected by the lasers. The occupation of orthogonal 
superpositions which can be excited by the lasers (often 
called bright states) decays due to the spontaneous emis- 
sion from the excited state. Many applications for optical 
systems followed, including laser cooling, lasing without 
inversion or coherent adiabatic passage [21 . 

As for some remarkable quantum optics phenomena, 
these applications have found their translation to meso- 
scopic electronic circuits There, transitions are medi- 
ated by coherent tunneling or by time dependent electro- 
magnetic fields. In quantum dot arrays, dark superpo- 
sitions are essential to proposals of coherent state trans- 
fer [llpr current switching by coherent population trap- 
ping ^. However, their detection in transport has been 
elusive for involving complicated structures 0, 0] or many 
level configurations Q. Importantly for this work, sin- 
gle electron spin resonance (ESR) has been achieved in 
quantum dot spin qubits driven by ac and dc magnetic 
fields [13] ■ Here we propose that the same setup can 
be used to detect a two electron dark state. 

Of great interest is the generation of entangled states. 
Entangling two qubits enables quantum information pro- 
tocols. In the solid state, entanglement has been demon- 
strated by the violation of Bell inequalities in supercon- 
ductor flux qubits [lH while quantum dot spin qubits 
have succeeded in the performance of essential ingredi- 
ents as universal quantum gates [13, 12 1 or single shot 
readout [l^]- Solid state entanglers have been proposed 
based on Cooper pair splitting [iJl, electron- hole exci- 



tations 15|, single electron emitters 16|, parity detec- 
tion [13] or purification protocols [l8|. Quite counter- 
intuitively, entanglement of macroscomc atomic ensem- 
bles can be generated by dissipation [l9| . A proper engi- 
neering of the environment drives an open quantum sys- 
tem to a pure steady state [13], similar to a dark state. 
In solid state qubits, the environment can be engi- 



neered by voltages and electromagnetic fields. Here we 
propose how to generate a maximally entangled dark 
state of two spatially separated electrons in a double 
quantum dot tunnel coupled in series to fermionic source 
and drain leads. In the Coulomb blockade regime with 
strong Coulomb interactions, electrons are transferred 
one by one through the system, which can be tuned to 
contain up to two conduction electrons, cf. Fig. [T] We 
consider a configuration described by the charge occupa- 
tion states (7VL,A^fl)=(0,l),(l,l),(0,2). In the presence 
of inhomogeneous in-plane magnetic fields, odd-parity 
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FIG. 1: Dynamical generation of entanglement, (a) In a spin 
blockade double quantum dot, charge flows through states of 
two electrons with opposite spins, {Itii); liit)}, coupled by 
intradot tunneling, r, to the doubly occupied singlet \Sr), 
from where an electron tunnels to the drain lead, (b) The 
occupation of the even-parity subspace blocks the current, 
forming a mixed state, psB, with statistical weights, Pa- Co- 
herent spin rotations by an ac magnetic field open the system 
to transport again. For a particular frequency, Q, an even- 
parity eigenstate of the magnetic field becomes a dark state, 
I'tpd)- (c) The system evolves towards a pure steady state com- 
posed of Bell superpositions. Coherence of the steady state 
is apparent in the undamped oscillations of the off-diagonal 
elements. 
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states, {It,!), lilt)}: support a current through the drain 
dot singlet |S'_R) = |0,ti) from which an electron is trans- 
ferred to the drain lead. However, the occupation of the 
even parity subspace, {|t,t), suppresses the cur- 

rent for a forward applied bias (from left to right) due to 
Pauli exclusion principle, what is known as spin block- 
ade Thus, the reduced density matrix of the double 
quantum dot, /5, evolves toward a mixed steady state in 
the even parity subspace: pss=Ptlt,t)(t,tl+Pili,l)U,ll, 
with statistical weights Single electron spin rotations 
by an ac magnetic field lift spin blockade and lead to a 
resonant current when its frequency matches one of the 
individual Zeeman splittings, ?iw=Ai,A/{ [13, 22. [23]. 
We find that for a particular frequency, superpositions of 
even-parity Bell states. 
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(|t,t)±li,i)), 



(1) 



decouple both from the driving field and transport. The 
system is thus dynamically driven to a pure dark state, 
Pst~\'4'd){'>pd\, see Eqs. (|4]) and ^ below. Any other su- 
perposition decays by the combined effect of the ac field 
and tunneling into the reservoirs. 

We investigate the fingerprints that such states leave 
in transport spectroscopy experiments as being carried 
out nowadays [lCll22l425|. Collective rotations of the two 
electron spins lead to a resonant current showing a sharp 
dip which cannot be explained by the individual electron 
spin dynamics. This effect will manifest clearly in the 
current at the border of the spin blockade window where 
the (1,1)— >-(0,2) tunneling transition is resonant. 

Model — Our system is modeled by a two sites An- 
derson Hamiltonian, i?=i?DQD+^fioads+.ffcoupi, describ- 
ing the double quantum dot, the electronic source (S) 
and drain {D) reservoirs, H\cads='}Zika'^ikd\k„dika, and 
their tunneling couphngs, .ffcoupi=I];fecr(Ai'^Lo-Cio-+H.c.), 
where the fermionic operators c^o- and dika annihilate an 
electron with spin a in dot iG{i, R] and lead l€{S^ I?}, 
respectively, and A; is the coupling strength. The dou- 
ble quantum dot term is given by -/?dqd — -^o~t~-^r5 ac- 
counting for the bare energy of the double quantum dot 
(including interactions), Hq^ and the coherent interdot 
tunneling, Hr=— X]o-('''Cl(tCrct+H.c.). We consider a sin- 
gle discrete level of energy £i in each dot which can be 
occupied by two electrons forming a spin singlet. Ex- 
cited states are assumed to be far off in energy so their 
contribution — forming on-site triplets that would lift 
spin blockade — can be neglected here. Thus we have 
^0 = SiCT e^cJcrCio- + Y.i Uifii^fiii + ULRULfiR, where Ui 
and Ulr describe on-site and interdot Coulomb repul- 
sion, and fiia and fii are the spin resolved and total num- 
ber operators, respectively. The chemical potentials of 
the leads, /x^, are such that only two electrons are allowed 
in the system: ei<fii—UhR<£i+Ui and Hi<ei-\-2ULR- If 
SR<fJ'D^ the right quantum dot stays occupied. Consid- 



ering spin, only seven states with charge distributions 
described by (0,1), (1,1) and (0,2) are relevant. 

Coherent interdot tunneling will be resonant when the 
states (1,1) and (0,2) are degenerate. We parametrize 
their detuning by 6e=eL—£R+UR—ULR- The spin block- 
ade window is hence defined by the region 6e > 0. 

For spin manipulation, we include magnetic fields in 
an ESR configuration [3, [3] ■ It consists of in-plane mag- 
netic fields with a dc component that creates a Zeeman 
splitting Af in each dot, and an ac component perpen- 
dicular to it whose frequency is close to the resonance 
conditions fku^A^ . Since magnetic fields are experimen- 
tally hard to localize, alternatives have been introduced 
using more tunable gate voltages. Effective ac magnetic 
fields are obtained by coupling the spin to a pulsed elec- 
tric signal mediated by either spin orbit or h ype rfine 
interaction 23-25[, slanting dc magnetic fields 22. 27|. 
or by spin-phonon coupling in the presence of mechani- 
cal vibrations of the dot [28| . For theoretical generality, 
we consider the purely magnetic field case given by the 
Hamiltonian term: HB(t)= X]i(^f cos wt, 0, Af )-Si, with 
M=9jifJ'BBji and spin operators Si=i J2aa' Ao^oa'ha'- 
The Lande factors, g^i, are in general inhomogeneous, [Ib 
is the Bohr magneton, B^i are the magnetic field j com- 
ponent in dot i, and & are the Pauli matrices. Different 
Zeeman splittings are required, Af T^Ap, otherwise the 
ESR field has no visible effect [26|. The different com- 
ponents A^ can be tuned either by the control of the 
Lande factor [2^ or by applying inhomogeneous magnetic 



fields: for instance by means of a micromagnet [22 



or pulses of different amplitude in electric dipole spin 
resonance experiments. Such experimental ability allows 
us to consider the simplest required configuration, with 
A|=A|j=7i3a; and A|=aA|j=7i?2. To ease the nota- 
tion, we include the asymmetries in the parameter a, as- 
suming that the magnetic fields are homogeneous. 

The dynamical evolution of the system is described by 
a markovian quantum master equation for the reduced 
density matrix, p, obtained by tracing the reservoir de- 
grees of freedom out: 



C 



VP, 



(2) 



where the commutator accounts for the coherent dynam- 
ics inside the double quantum dot. The Liouvillian su- 
peroperator ^T=^a=±i^f describes tunneling events 
to (-I-) or from (— ) lead I and includes decoherence for 
the finite lifetime of states coupled to the leads [2^ . 
Spin decoherence is assumed to be of a longer time scale. 
The drain current is then given by I—qti[{C]^ — C~j^)p\. In 
the high bias regime, Cg—Cj^—Q. The remaining terms 
are fully described by the tunneling rates Ti = '^\Ki\^vi, 
where vi is the density of states in the leads, for pro- 
cesses carrying an electron from the source to the dots 
and from the dots to the drain. Processes taking elec- 
trons in the backwards direction are negligible, though 
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FIG. 2: Current as a function of the driving frequency, uj. 
Collective spin rotation results in a nontrivial lineshape not 
centered at the ESR conditions (vertical dashed lines) and 
a vanishing current (dark state) at ?if2=(Af,-|-A|{)/2. Cases 
for different driving amplitudes, B^, are vertically offset. Pa- 
rameters: ftr=10"^meV, r=20r, 7=0.04meV/T, B2=1T and 
a=0.8. Interdot tunneling is resonant: fc=0. The inset shows 
the shift of the antiresonance frequency, Um, with the ampli- 
tude of the ac field. The solid line is a quadratic fit. 



they are included in the numerical calculations. 

Results — Naively, one would expect that the current 
spectroscopy results from the sum of two Lorentzian res- 
onance peaks centered at fta;=A|^, A|j — when the rota- 
tion of each electron spin would lift spin blockade — and 
is zero elsewhere. Indeed, this picture agrees with exper- 
imental observations 2^, 24, 25|. However, due to the 
interplay of coherent spin rotation and interdot tunnel- 
ing, resonant features of a non trivial lineshape appear, 
as shown in Fig. O On one hand, peaks are not nec- 
essarily centered at the individual spin resonances but 
at points of maximal hybridization which depend on t, 
and 5£. On the other hand, and most importantly, 
current vanishes for a frequency fta;=ftri=(A£-f A|j)/2. 
Such anti-resonant behaviour is a clear signature of dark 
states. We note that similar current characteristics are 
predicted in discrete lattice models of transport 13(| to 
which our system can be mapped. 

Analytical understanding of the dynamics at that point 
can be obtained within a rotating frame approximation, 
when Af <cA|'-^?za;. Thus, neglecting the contribution of 
counter-rotating terms, we get a time independent mag- 
netic field Hamiltonian: 



(3) 



We can easily verify that, at w=ri, |?/'.-) = (|tt)-|ll))/\/2 
is a zero eigenvalue eigenstate of HB,wNA+Hr- Therefore 
it is decoupled both from the external magnetic field and 
from transport, due to parity and spin blockade. Elec- 
trons occupying any other state or superposition will flow 
to the drain lead and be replaced. As a consequence, 
the system will dynamically evolve towards a stationary 




FIG. 3: A pure coherent steady state, (a) Condition for 
a pure state as a function of the driving frequency. Solid 
lines represent the diagonal elements of the density matrix. 
Green dots are for off-diagonal terms. At the antiresonance 
frequency, hfl, a pure state is formed, p^~p- Here, Bx~Bz/5. 
(b) State tomography at uj—Q. Only two-electron states are 
represented for clarity. Oscillatory off-diagonal elements are 
averaged out but will affect the variance of the measurement. 



solution given by the Bell state The density ma- 

trix will be then described by the pure state 
fulfilling £rwa|''/'-) (V'- 1=0, with the Liouvillian super- 
operator CmvAO=-ifi^^[HBQu+HB,mvA,0]+Crd. In 
this sense, opening the system to transport drives it to a 
maximally entangled dark state, for any initial state. 

In the laboratory frame, the steady state describes 
Rabi oscillations of the two even-parity Bell states, |V'±)j 



with frequency f7=(A|,-f A|j)/(2fi.): 

\4'stit))—i sinilt|^+)-(- cos iltlip^) 



(4) 



within the rotating wave approximation. Hence we 
can approximate the stationary density matrix by 
\ipst{t)){'(l^st{t)\, explicitly: 



/5.tW«^[lt,t>(t,tl + IU>(U|-(e 



'''"it,t>a,;i+H.c.)]. 

(5) 

We note that this solution is exact for a circularly polar- 
ized ac magnetic field but, for experimental feasability, 
we restrict our analysis to linear polarization. As known 
since the work of Bloch and Siegert [sH , counter-rotating 
terms of linearly polarized electromagnetic fields induce 
a shift of the resonance condition, see Fig. [5] 
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FIG. 4: Transport spectroscopy for different dc magnetic 
fields, witfi Bx=lT- (a) For finite detuning {Se=0.2 meV), a 
double peak structure is consistent with previous experimen- 
tal observations (see text), (b) Resonant interdot tunneling 
((5£=0) enhances the visibility of the collective features. The 
dark antiresonance is robust upon increasing the magnetic 
field. An offset is introduced for clarity. 



A numerical confirmation of our above results is shown 
next. As every pure state, the dark state must be de- 
scribed by an idempotent density operator: p^—p. Oth- 
erwise it is a mixed state. In Fig.[3Ja) we verify this prop- 
erty for the stationary solution of the full time-dependent 
master equation ([2]). There, the steady state is mixed, as 
expected for a transport configuration. However, at the 
antiresonance frequency it is dynamically purified. 

A detection scheme for state coherence is quantum 
state tomography. Single-shot readout of the two elec- 
tron spins can extract the required correlations In 
our case, the explicit time dependence of the dark state 
non diagonal elements, cf. Eq. ([S]), will be averaged out 
in any measurement, see Fig. ^h). The result would 
be then similar to that of a mixed spin blockade state, 
PSB=Scr?'o-|o',(T)(cr,(T|. However, entanglement of the 
dark state will be explicit in the variance of the mea- 
surement, see Fig. [2Ub). 

Experimental discussion — No signature of such a 
dark Bel l state has been found in the experiments up 
to date 



22|-|25l|. Our results indicate that the role of 



coherent tunneling must be emphasized. Usually, config- 
urations where interdot tunneling is not resonant show 
well defined double peak resonances centered at the in- 
dividual ESR conditions, fia;=A|, A|j |22| - |25{ . Using re- 



alistic parameters, our model reproduces that behaviour 
for fe>0, cf. Fig. mja). In that case, the effect of coher- 
ent tunneling is reduced and only a tiny antiresonance 
appears. 

We propose that a clearer evidence of the dark Bell 
state will appear in the resonant tunneling case, ^£=0, 
due to the enhanced interplay of coherent tunneling and 
collective spin resonance, see Fig. Hfb) . For low , two 
peaks appear around the antiresonance condition, whose 



positions weakly depend on the magnetic field and are not 
centered at huj=AL, Ar. There, the sharp antiresonance 
structure is dominant, with a vanishing current. As an- 
alyzed in Ref. 0, a finite current at the dark state con- 
dition can be used to estimate the effect of other sources 
of decoherence different from tunneling. Upon increasing 
the dc magnetic field, the current develops humps that 
follow the individual ESR conditions and would be even- 
tually resolved as separate peaks for larger B^. However, 
the central structure is robust to decoherence and con- 
forms an unambiguous signature of the dark state, even 
if the antiresonance becomes faint. 

Conclusions — We have predicted a transport induced 
maximally entangled state of two spatially separated 
electrons in a driven double quantum dot. Interplay of 
coherent interdot tunneling and collective electron spin 
resonance is essential for leaving clear signatures of the 
dark Bell state in transport spectroscopy experiments. 
The entangled state is decoupled from its environment 
by parity symmetry and Pauli exclusion principle. Thus 
it is not affected by decoherence due to coupling the sys- 
tem to leads. Other sources of decoherence can be probed 
and will motivate further investigation. We propose an 
experimental configuration which is within reach where 
resonant interdot tunneling (at the border of the spin 
blockade window) enhances the visibility of the entangled 
dark state features. Our work introduces a mechanism to 
produce Bell states in open systems for any initial condi- 
tion. It will allow for investigations of non local quantum 
correlations of two electrons stored in solid state qubits. 
Our simple configuration constitutes an ideal candidate 
for the detection of transport dark states. 
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